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1. Introduction 

The Lentz-Ising model of ferromagnetism has been thoroughly studied since its 
conception in the 1920's[Lcn20j. It was solved in the 1-dimensional case by Ising 
himself in 1925[Isi25 and in the 2-dimensional case without an external field by 
Onsager in 1944 |Ons4 4 . For an introduction to the model, see Cipra |Cip87| . 

The partition function on a graph G on n vertices and to edges is defined as 

(1) ZK/: ,-.//) ^>,...,-y. 

Here ay counts the number of induced subgraphs of G with (n — j) /2 vertices and 
(to — i)/2 edges in the boundary. We refer to the index i as the energy and the 
index j as the magnetization. 

The traditional partition function studied in statistical physics is then obtained 
by evaluating it at a certain point 

(2) Z(G;e K ,e H ) 

where K — —J/ksT, H = —h/ksT and J and h are parameters describing the 
interaction through edges and with an external magnetic field respectively, T is the 
temperature and ks the Boltzmann constant. 

The main goal in the study of the Ising model on a graph G, or some family 
of graphs, is usually to study the model in the vicinity of a critical temperature, 
denoted T c , where the model undergoes a phase transition and there determine the 
behaviour of various critical properties. 

2. Definitions and notation 

Let G = (V, E) be a graph with vertex set V with \V\ = n vertices and edge set 
E with \E\ = to edges. Let a state a be a function from the vertices V to the set 
{±1} and let be the set of all states. We can then define the energy of the graph 
G = (V, E) in state a to be 

(3) S(G,a)= °(u)a(v) 

uv&E 

the magnetization to be 

(4) M(G,a)=Y J <v) 

and formulate a generating function that counts them all as 

(5) Z(G; x,y) = J2 x £ ^y M ^ = £ «,,<'V 

where the last equality defines the coefficients ay — aij (G) . We often drop G when 
we can deduce the graph from the context. 
We will also need these definitions later: 
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Definition 1 (T-join). A T-join (T,A) in a graph G = (V,E) is a subset T C V 
of vertices and a subset A C E of edges such that each vertex in T is incident with 
an odd number of edges in A and each vertex in V \ T is incident with an even 
number of edges from A. 

Observe that the cardinality of T has to be even since we can not have a subgraph 
with an odd number of vertices of odd degree. 

Definition 2 (Cut). A cut [S, S] in a graph G — (V, E) is a subset of edges, induced 
by a partition S U S = V , that have one endpoint in S and the other in S. Let 
| [S, S] | be the number of edges in the cut. 

Definition 3 (Locally vertex transitive). We say that a sequence of graphs {Gi}°l 1 
are locally vertex transitive if for each R there exists an N such that all balls of 
radius R, in the graph metric, around each vertex in each graph in the subsequence 
{Gi}°^ N are isomorphic. 

3. Series expansion 
We can try to make our problem simpler by setting y — 1 and get 

(6) Z(x, 1) = Z(x) = OijafV = Y a ^ 

ij i 

where once again the last equality is the definition of the <ij coefficients. What do 
these coefficients count? Since the state a partitions the vertex set V in two parts 
and the £(G,a) counts edges with one vertex in one part and the other vertex in 
the other part as negative and the other edges as positive we get that en is twice 
the number of cuts of size m f 1 (we count each cut twice since we can interchange 
the partitions). This has a natural reformulation using even subgraphs, namely: 

Theorem 4 (van der Waarden). Let at be the number of cuts of size and let 
bi be the number of even subgraphs with i edges, then: 

(7) j2 a ^ lK = 2coshm K J2 btt&nhl K 

i i 

Proof. The first sum in Q is the moment generating function for the sequence a^. 
The k th moment of can be written as 

(8) Mfc = y ^ k = Yi ( Yi ^H^H ) 

i <tGO \uvEE / 

We now expand the multinomial {^2 uveE o~(u)o~(v)j where each term can be seen 
as an choice of k out of m edges (not necessarily distinct). Now observe that if we 
have chosen an even number of edges incident with a vertex, v say, we will have an 
even number of a(v)'s in the product so they contribute +1. If we have chosen an 
odd number of vertices we can find a smallest (in some arbitrarily order) such odd 
vertex v and we see that if we change the state a to the state a' with a(v) = —a'(v) 
and all other values equal, we will get a bijection between states witch contribute 
+ 1 and —1 and with at least one vertex of odd degree. Our conclusion is that 
we only count the choices where we have an even degree at each vertex. We will 
however count subgraphs where we have the opportunity to choose each edge a 
multiple number of times. If we reduce the multiple edges modulo 2 we get a 
simple subgraph of even degree. The "surviving" edges are the ones that where 
chosen an odd number of times so an even number of those "odd" edges have to be 
incident at each vertex. 
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If we now change our view and instead of adding up the k moments, change the 
order of summation, and add up along the index of the number of "surviving" odd 
edges i we get a simple connection between the simple subgraphs and subgraphs 
with multiple edges. We can construct an even multiedge subgraph by first select 
a simple subgraph with even degree at each vertex and then multiply each edge an 
odd number of times and then select a number of edges not in the even subgraph 
to multiply an even number of times. So if we first choose an even subgraph with i 
edges and multiply each edge an odd number of times we get the generating func- 
tion bi sinh 1 K, and then choose a number of edges outside the even subgraph and 
multiply these an even number of times, we get the generating function cosh m ~ l K, 
and we end up with an multiedge subgraph with an even number of edges incident 
to each vertex. If we now sum over all i we get 

(9) ^ h sinh J K cosh" 1 " 1 K = cosh" 1 K^b, tanh 1 K 

i i 

and we see each of these graphs twice and therefor this gives Q . □ 

To formulate the full two variable connection we need T-joins instead of even 
degree subgraphs and also consider the size of the sets in the vertex partition 
induced by the state. Also note that in this theorem we have slightly changed the 
meaning of and bij to make the proof using standard graph theoretic notation. 
The following theorem can be found in e.g. |Big77| : 

Theorem 5. Let G = G(V,E) be a graph, the number of cuts [S,S] with 
\[S,S]\ — i and \S\ — j. Let bij be the number of T-joins {T,A) with \A\ = i and 
\T\ = j . Then 

MV = 2~ |l/| ]T %(i - x y(i + .t)I £ i-*(i - y y(i + y)\ v \-i 

ij ij 

Proof. Fix a subset T C V of vertices and a subset A C E of edges from the graph 
G = G(V, E). Let S C V be another subset of vertices and [S, S] be the cut defined 
by the edges from S to S = V \ S. Let the weight of the vertices in T n S be — y, 
the weight of the vertices in T n S be y, the weight of the edges from A that lies in 
the cut [S, S] be —x and the rest of the edges from A have weight x. Let the total 
weight of (T, A) with respect to the cut [S, S] be the product of the weights of the 
edges and vertices in (T, A) . We say that the weight is positive if the coefficient in 
front of a;'" 4 'y' T ' is positive and negative otherwise. By magnitude we denote the 
weight without the sign. 

(T, A) can fail to be a T-join in basically three ways. First the cardinality of 
T can be odd, secondly there can exist a smallest vertex (in an arbitrary order of 
the vertices) v that is incident with an odd number of edges from A and does not 
belong to T, and finally there can exist a smallest vertex v that is incident with an 
even number of edges from A and belongs to T. 

In the first case we have two cuts [S, S] and [S, S] in which the magnitude of the 
weight will be the same but with opposite sign. 

In the two latter cases we have a bijection between cuts with v E S and v £ S 
(we simply move the vertex v between S and S) that once again give the same 
magnitude and different signs of the weight. If we sum over all cuts the total 
contribution of such a choice of (T, A) will cancel. 

If (T, A) indeed is a T-join the weight will always be positive since we either have 
an even number of vertices in T n S and an even number of edges crossing the cut 
or an odd number of vertices in T D S and an odd number of vertices crossing the 
cut. All in all we end up with an even number of minus signs and thus a positive 
weight. 
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If we now sum over all choices (T, A) and S we will count each T-join 2' y l times. 
If we rearrange our summation (i.e. we first choose a cut and then go through all 
choices of T and A) we get the theorem. □ 

To get Theorem |4] we have to shift the indices and substitute x for e~ 2K and y 
for and finally double all values since we count all cuts twice in ([7]). 

(10) 2e mK J2 a>ije~ 2iK & = [since 0° = 1] = 2a^ m - 2l)K = 

2e mK 2- m J2 M 1 - e- 2K ) l (l + e - 2K ) m ~n n = 

2 cosh™ K h tanh 4 K 

i 

where bi denotes the number of T-joins with i edges. Now, since counts the 
number of cuts with i edges in Theorem|H] 2a.;e ( - m_21 ^ corresponds to the aie lK in 
Theorem [4] via reindexing and the fact that each cut is counted twice in Theorem 
51 

Since we have a symmetry between T-joins and cuts we have the following corol- 
lary: 

Corollary 6. With the same notation as in theorem^we have 



Proof. If we choose a T-join instead of a cut the weight of (T, A) will always be 
positive if and only if (T, A) is a cut. In other cases the contributions once again 
cancel out. □ 



3.1. The thermodynamic limit. In Physics we are interested in the so called 
thermodynamic limit of a sequence of graphs {Gj}^ l . This is defined as 

(11) f{x)= \\m oW ^ m \ogZ{G 3 -x,l) 

when it exists. An example of such a family is {C n x C„}^L 3 . If the graph family 
is locally vertex transitive its easy to see that the number of connected T-joins 
of fixed size will grow proportionally to |V^(G)|. From that follows that the total 
number of T-joins of a fixed size will grow as a polynomial with degree equal to the 
maximal number of connected components in the T-joins of that size and thus will 
the thermodynamic limit exist. 

If we change notation so that our index set instead is the number of vertices n 
in our graph sequence and use bj(n) to denote the number of T-joins with j edges 
we see that the thermodynamic limit is 

(12) lim ^ log 2 cosh" 1 x bj (n) tanh J x — log 2 + d cosh x + bj tanh- 7 x 

j 3 

which defines a new set of bj'.s that happens to be rational numbers and d — m/n, 
the average degree. 
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3.2. Taylor series. It can be of interest to plot these functions, or at least some 
approximation of them. Since we have a phase transition in all interesting cases its 
hard to find one function that works for the entire interval. In this section we will 
instead develop two Taylor approximations, one for the high-temperature case and 
one for the low-temperature case. The first one is simple, take 



(13) 



i(K) =V a t K l =log2 + rflogcoshif+ lim ± log V b 2 (n) tanh 4 K 



and Taylor expand around K = 0. If you want the function K(u) you can invert 
u(K) as a Taylor scries. Its often better to plot a Pade-approximation of u(K). 
The second one needs a little more work: 



(14) " 



dK 



log 2 + d log cosh K + lim - log b t (n) tanh* K 



x = tanh K 



dx - 



(1-x 2 ) — lim Mog YhMx 1 = 



dx + (l- x 2 )— Jim i log £ ^ a m _ 2i (n)(l - x)*(l + or)' 



g?x + (1 — a; 2 )— lim 

dx n— >oo 



-log 2 + dlog(l + x) + Mogexp ( n^ai (j 



x 



x 

1 - X 



X 



(x+1) 



1 -x 
1 + x 



„ \ ^ / 1 — X 
d- 2 > ai 



d 



- 2 Em t 



1 - tanh K 



tanh jY 



By not doing the last substitution x = tanh K its easy to invert this function and 
plot it in the low temperature (high energy) portion of the scale. 



4. TWO OTHER SERIES EXPANSIONS 

From the basic thermodynamic limit one can construct two other sequences that 
are of a more combinatorial flavour. These are sequences with integer coefficients. 
This is because they, instead of weighting the counts, make an explicit order in 
which you have to choose things and thus avoid dividing with large factorials. The 
correspondence with the thermodynamic limit series are: 

(15) 5>(ny = exp (nJ2hA = ]J i^)"' = II (tZW^Y 
i>a \ i>i J »>i v ' 7 i>i V1 Pi x ' 

The coefficients are connected trough these simple triangular linear equations: 

(16) ib i = J2t/3t 

t\i 

(17) ib t = ]T rfi 

rs—i 

(18) ^/3 ^ = ]T Mr)t& 
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4.1. What they are counting. The new series can be seen as placing connected 
subgraphs at each vertex in some order. First we place all subgraphs (including the 
empty one) rooted at vertex one, then the subgraphs rooted at vertex two and so 
on in such a way that we never form a new connected component. In this way we 
avoid symmetries and we get a integer sequence. We can still calculate, in principle, 
the different /3j and (3i for locally vertex transitive graphs. 

5. Plots 

We shall now compare these Taylor expansions of the series with some sampled 
data to compare how well behaved the series are around the critical energy. As we 
shall see, the series are rather far from what can be expected to be the truth. 




0.2 0.4 0.6 0.8 1 



Figure 1. K(u) for the three-dimensional simple cubic lattice to- 
gether with sampled data for the cubes of linear order 32, 64 and 
128. 

5.1. Simple cubic lattice. We start of with the simple cubic lattice in three, 
four and five dimensions. As can be expected, the longest series expansion is for 
the three-dimensional simple cubic lattice. The four- and five-dimensional simple 
cubic lattices have much shorter high- and low-temperature expansions. Figure 
[T] to [3] shows K(u) curves for the three lattices respectively. The pictures are of 
diagonal Pade-approximants of the Taylor expansions. The Pade-approximant for 
the function of the three-dimensional simple cubic lattice is fairly accurate to about 
K w .22. For the four- and five-dimensional lattices the accuracy is a lot lower as 
can be seen from the pictures. 

6. Guessing the radius of convergence 

The interesting part is of course to try to find the radius of convergence for 
the series of the function u(K) since that would give the critical temperature K c . 
A fairly good guess for K c for the simple cubic lattice is that it is approximately 
0.2216546 (see e.g. [HRL + ]). A simple observation is that the radius of convergence 

. . . . J I 

must be smaller than any l//3j since otherwise the infinite product (|15| will not 
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Figure 2. K(u) for the four-dimensional simple cubic lattice to- 
gether with sampled data for the cube of linear order 4, 6, 8, 12, 
16, 32, 48 and 64. 

0.5 r 



0.4 - 



0.3 - 
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Figure 3. K(u) for the five-dimensional simple cubic lattice to- 
gether with sampled data for the cube of linear order 4, 6, 8, 12, 
16 and 32. 

converge. In Table [T] the row labelled min is the minimum over all known values of 
(3i for the simple cubic lattice. In trying to extrapolate this value we have used a 

^ l/i 

polynomial of degree 3 and a min square approximation to the data set (l/i, ) 
and looked at the constant term in the resulting polynomial. This gives the result in 
the row labelled /(0). Since the coefficients is some form of combinatorial quantity 
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A 


h 




a~ 


4 


« 4 


min 


0.291989 


0.291989 


0.618531 


0.616299 


0.618531 


0.616307 


/(O) 


0.227727 


0.227839 


0.54278 


0.543852 


0.542846 


0.539196 




0.221418 


0.221451 


0.523726 


0.520324 


0.523747 


0.542904 




0.258385 


0.258429 


0.578131 


0.565115 


0.578137 


0.56453 



Table 1. Different ways to try to guess the radius of convergence 



for the simple cubic lattice. 



it is not unthinkable that they grow exponentially. We have thus also tried to 
fit the data (i,(3i) to the functions a l i^^i and a lr y and calculated 1/a that gives 
us the guessed radius of convergence. This is the last two lines of Table [T] The 
second column is the same thing for the coefficients of bi. When we come to the low 
temperature part we get into some trouble since the coefficients of the sequences 
Oi and cti have both positive and negative signs. If we try to use the same analysis 
as for the (all positive) bi and A, we get very erratic numbers, so instead we split 
the sequences in a positive and a negative part and do the analysis separately. As 
observed by others, the low temperature series seems to have a complex root that 
is closer to the origin than the physically important real root. Fortunately does 
the extrapolation with a l i°^ for the bi and A series give a decent idea of what the 
critical temperature may be. 
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Appendix A. Tables 

We have collected the various series we have found in this appendix. Some of 
these are old and thus not especially long and some are from newer calculations 
and longer. 



n 


a n 


a n 


a n 


6 


1 


1 


1 


10 


3 


3 


3 


12 


-4 


-4 


-7/2 


14 


15 


15 


15 


16 


-33 


-33 


-33 


18 


104 


104 


313/3 


20 


-282 


-285 


-561/2 


22 


849 


849 


849 


24 


-2460 


-2470 


-9847/4 


26 


7485 


7485 


7485 


28 


-22542 


-22647 


-45069/2 


30 


69392 


69384 


346966/5 


32 


-213738 


-214299 


-427509/2 


34 


666750 


666750 


666750 


36 


-2086785 


-2092121 


-12520405/6 


38 


6583341 


6583341 


6583341 


40 


-20852223 


-20892996 


-83409453/4 


42 


66425750 


66424630 


464980286/7 


44 


-212410377 


-212770353 


-424819905/2 


46 


682202205 


682202205 


682202205 


48 


-2198562644 


-2201602421 


-17588511087/8 


50 


7110521070 


7110521022 


35552605353/5 


52 


-23065955826 


-23093964696 


-46131904167/2 


54 


75045653088 


75045278168 


675410878105/9 


56 


-244806881325 


-245063348553 


-979227570369/4 


58 


800606679471 


800606679471 


800606679471 


60 


-2624325216574 


-2626724535242 


-13121625909861/5 


62 


8621219166681 


8621219166681 


8621219166681 


64 


-28379404026078 


-28402366460136 


-113517616531821/4 



Table 2. Simple cubic lattice 
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SERIES EXPANSION FOR THE DENSITY OF STATES OF THE ISING AND POTTS MODEI1S 



n 


<4 


K 


a 5 „ 




4 





6 





10 


6 





76 





180 


8 


1 


1371 





5025 


10 





30152 


1 


178696 


18 







5 




20 


28 




-11/2 




22 


-64 




45 




24 


127/3 




-95 




26 


228 




50 




28 


-834 




5 




30 


1116 




1471/3 




32 


5911/4 




-1685 




34 


-10404 




1885 




36 


21460 




-775/2 




38 


-1956 




5445 




40 


-595179/5 




-112951/4 




42 


1076092/3 




49545 




44 


-344316 




-62795/2 




46 


-1132588 




71442 




48 


10842287/2 








50 


-9187444 








52 


-5820150 








54 


73867260 








56 


-1294335811/7 








58 


95069292 








60 


2609680726/3 








62 


-3217644924 








Table 4. 


Simple cubic lattice in dimension 4 (x^) and 
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DANIEL ANDREN 



n 


a n 


K 


4 





6 


6 


2 


44 


7 





36 


8 





384 


9 





688 


10 


6 


4572 


11 


6 


11148 


12 


-14 


66158 


13 





190662 


14 


30 


1051668 


15 


60 


3471452 


16 


-108 


17917704 


17 


-144 


65438160 


18 


1118/3 


971627006/3 


19 


498 


1265584728 


20 


-714 


30625029636/5 


21 


-2366 


25078631014 


22 


3270 


-3816568476 


23 


7704 




24 


-8106 




25 


-27372 




26 


19842 




27 


114342 




28 


-68892 




29 


-377376 




30 


643952/5 




31 


1431726 




32 


-137718 




33 


-5365756 




34 


-33330 




35 


18644574 




36 


-970922/3 




37 


-69012330 




38 


-32516754 




39 


249820316 




40 


162829320 




41 


-869879742 




42 


-5660822830/7 




43 


3155460756 





Table 5. Simple cubic lattice, 3-states Potts model 



